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Abstract 

This paper presents the Thomas-Fermi approach generalized to consider the 
particle correlations in many-body systems with non-Coulomb interaction poten- 
tials. The key points of the generalization consist in using integral formulation and 
extracting the pair distribution function. The latter has been found to obey the in- 
tegral equation which, in the classical limit, is reduced to the well-known equation of 
Bogoliubov. So, the approach presented gives the way of taking into account inter- 



actions in the quantum many-body systems beyond the weak coupling perturbation 
theory. 

It is known that Thomas-Fermi approach [jlj together with its particular 
variant for the classical many-body systems, the Debye-Hueckel method |2j], 
is able to provide satisfactory results while investigating the correlation 
phenomena in the media with Coulomb interparticle potential. The advan- 
tage of this approach consists in the extreme transparency of understanding 
the physics of the particle correlations. But up to now the applications of 
Thomas-Fermi method have been limited to the situation of the Coulomb 
type interactions. This paper presents the generalization of the approach 
to consider the many-body systems with constituents interacting by means 
of an arbitrary relevant potential. 

To derive the generalization, let us take the system of some particles 
uniformly distributed with density n over a region of volume V. Let the 
potential of the particle interaction be $(r) . Assume that a point object 
is placed into this region and interacts with the surroundings via the same 
potential. Further we shall treat the point of the object location as the 
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origin of the coordinates. We are interested in the equilibrium structure 
produced around the origin point. This structure is specified by quantity 
n str (r;9,n) being the density of the particles located in the vicinity of 
point f , r =| f |. Here 9 denotes the system temperature. Below we shall 
use n str {r) instead of n sir (r;#, n) to simplify formulae. According to the 
first "item" of Thomas- Fermi approach, n str {r) can be calculated with the 
following condition: 



where Pid(p) is the chemical potential of the ideal gas of the considered par- 
ticles at density p, and U (r) is the energy of the interaction of the particle 
being at point f with the other particles of the system. The quantity in 
the left side of (1) can be named the total particle potential. To calculate 
n s tr(r) with (1), the connection of U(r) with n str (r) should be fixed. In the 
case of the systems with Coulomb interparticle potential, this is realized 
by employing the neutralizing background and Poisson equation [1,2]. But 
it is impossible for the case of non-Coulomb interactions. However, we can 




(i) 
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follow the other way based on the obvious integral relation 

U(r) = J v mf-y\)n str (y)dy. (2) 

At point ff ar which is far enough from the origin of the coordinates, the 
total particle potential then takes the value 

Hidiji) + n j v ${\r far - y\)dy (3) 

if, of course, we hold to the reasonable assumption that 

lim $(r) = . 



Note, that we consider rj ar = | rjar |<C R, where R is the minimal distance 
between the location point of the object and the system boundary. So, any 
boundary effects can be ignored in our investigation. Equating (1) and (3) 
we can find 

tHd(n) - Hid (n str (r)) = $(r) + j y (n str (y) - n)$(| r-y \)dy , (4) 
where the used relation 

lim Svmr fa r-y\)-^\r-y\))dy = (g) 

v^o J v <S>(\r-y\)dy 
is fulfilled for a lot of the known potentials, in particular, for the Coulomb 
potential as well as for the integrable ones. 



Now we should realize the second "item" of Thomas- Fermi method. It 
implies that using n str (r) : we are able to study the space correlations in 
the many-body system made of the particles uniformly distributed with 
density n. The most important quantity calculated in Thomas-Fermi ap- 
proach and related with the particle correlations, is usually thought to be 
the screening radius. Of course, now we can not exploit this quantity, for 
it is tightly connected with Coulomb interparticle potential and the neu- 
tralizing background. But there is one more correlation characteristic, no 
less important than the screening radius, that may be evaluated in the 
approach discussed. It is the pair distribution function g(r; 6, n) for which 
we have 



Using (4) and (6) we readily obtain the following integral equation for the 
pair distribution function: 



where g{r) = g(r; 9, n). To receive some view to what extent this equation 




(6) 



lkd{n) - Hid (ng(r)) = $(r) + n J y (g(y) - 1)$(| f-y \)dy , (7) 
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is adequate, let us explore its classical limit. In this case we have 

tHd(n) - Hid(ng(r)) = -9 In g(r). (8) 

Therefore, relation (7) can be rewritten as 

-9 In g{r) = $(r) + n J v {g{y) - 1)$(| f - y \)dy . (9) 

Equation (9) makes it possible to generate the expansion in powers of n 
for g(r; 9, n): 

g (r; 6, n) = g {r; 9) + n gi(r; 6) + . . . , 

where 

g (r; 9) = exp 

9i(r;0) = --expl ^1 j v ^{\r-y\) (exp( f±) - 1) dy . 

As it is seen, go(r; 9) is in full agreement with the known result of the 
calculations in Gibbs canonical ensemble @. But expression for gi(r;9) 
coincides with the true value 

$(r)' 



~9~ 



gi(r;9) = exp 
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J v (exp(- Hlf ~^ l) ) - 1) (exp(-^f) - l)dy 



only at high temperatures. So, equation (9) gives the valid second virial 
coefficient at low temperatures, and the correct second and third ones at 
high temperatures. Note, that quite reasonable integral equations for g(r) 
derived in the superposition approximation Hfl, are in error beyond the 
third virial coefficient [|3], [5|]. Thus, the use of integral formulation for 
Thomas-Fermi approach instead of the differential one and the orientation 
to the pair distribution function, enable us to generalize this method to 
study the case of the non-Coulomb interactions. But the investigation goal 
has not been reached yet. 

Indeed, equation (7) gives the possibility to consider particle correla- 
tions in, for example, nuclear matter, where the nucleon-nucleon potential 
is of Yukawa type exp(—y)/y, or in the electron gas. But it is not the 
case for the potentials which are not integrable due to their behaviour at 
small separations between interacting particles. For instance, if we take 
potential \jr m (m > 2) then the integrals in expressions (7) and (9) will 
not exist. This obstacle appears because particle correlations have been 
neglected while calculating U(r). In particular, in above mentioned reason- 
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ings, quantity U(rj ar ) is given with the expression 

U{r far ) = n J v $(\ f far -y\)dy 

that corresponds to the particle interactions taken in Hartree approxima- 
tion. The latter is well-known not to take into account any correlations. To 
derive a more correct relation for quantity U(rf ar ), note, that it is noth- 
ing else but the interaction energy of a particle with the other particles 
uniformly distributed around with density n. Keeping this in mind, we 
estimate U(rf ar ) [J6j as 

U(r f ar)=nJ v g(\ f far -y\) $(| f far -y\)dy. (10) 

But now the problem arises how U (r) can be specified, for we have fixed 
only its limit value at ff ar . A reasonable way to do this is to treat quantity 
g(r)$(r) as the effective interaction potential which should be substituted 
for $(r) in all the integrals in expressions (2)-(5). In this case we have 

U(r) = J y g{\ r~y\) $(| r-y\) n str {y)dy . (11) 
Further, (1) and (11) together with (6) yield the equation 

tMd(n)-fHd(ng(r)) = $(r)+n j y (g(y)-l) g{\ r-y\) $(| r-y\)dy . (12) 
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It is interesting to compare (12) with the well-known integral equations 
for the pair distribution function of the classical simple liquids. In such a 
situation expression (12) can be rewritten as 

- 9 In g{r) = $(r) + n J v {g{y) ~ 1) g{\ f- y |) $(| r - y\)dy. (13) 

This is very similar to the equation 

-6 In g{r) = $(r) + n J y {g{y) - 1) 5 (| f - y \)dy (14) 

derived in the superposition approximation by Bogoliubov [[5j . Here 

5 (I r- y |) = g{\ f- y |)$(| f - y |) + F $(t)^-dt. 

J \r-y\ dt 

Function 3>(t)^f is oscillating around zero, so its integral, apparently, does 
not make an essential contribution into the expression for $ (| r — y |) . So, 
equations (13) and (14) have to yield the similar qualitative results. Ac- 
cording to (14), we can expect that $ (r) is more accurate estimation of the 
effective interaction potential than g(r)<f>(r). The replacement of g(r)<fr(r) 
by $ (r) in expression (12), results in the following integral equation: 

Vid(ri) - md(ng(r)) = $(r) + n J y (g(y) - 1) 5 (| r-y\)dy, (15) 
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which gives Bogoliubov equation (14) in the classical limit. Thus, we have 
now derived integral equations (13) and (15) that allow us to explore par- 
ticle correlations even in the many-body systems with the interaction po- 
tentials behaving at small distances as \/r m (m > 2) . 

In conclusion, let us take notice of the most important points of the 
paper once more. The article presents the generalization of the Thomas- 
Fermi method to investigate particle correlations in the many-body sys- 
tems with non-Coulomb interactions of their constituents. The essentials 
of the generalization consist in using integral formulation and operating 
with the pair distribution function. The paper results in the integral equa- 
tions for g(r) which correspond to the different ways of considering particle 
interactions. The most accurate of them, (15), is reduced to Bogoliubov 
equation (14) in the classical limit. As to relation (7), it follows from (15) in 
approximation g(r) « 1 . At last, integral equation (12) can be used as the 
simplified variant of (15) being able to produce the same qualitative picture 
of particle correlations. Of course, it is also interesting to investigate other 
possible evaluations of U (r) and the corresponding integral equations for 

the pair distribution function. This investigation will be presented in the 
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forthcoming paper. 
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